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We present a complete theory of the spin torque phenomena in a ultrasmall nanomagnet coupled 
to non-collinear ferromagnetic electrodes through tunnelling junctions. This model system can 
be described by a simple microscopic model which captures many physical effects characteristic 
of spintronics: tunneling magneto resistance, intrinsic and transport induced magnetic relaxation, 
current induced magnetization reversal and spin accumulation. Treating on the same footing the 
magnetic and transport degrees of freedom, we arrive at a closed equation for the time evolution 
of the magnetization. This equation is very close to the Landau-Lifshitz-Gilbert equation used in 
spin valves structures. We discuss how the presence of the Coulomb blockade phenomena and the 
discretization of the one-body spectrum gives some additional features to the current induced spin 
torque. Depending on the regime, the dynamic induced by the coupling to electrode can be viewed 
either as a spin torque or as a relaxation process. In addition to the possibility of stabilizing uniform 
spin precession states, we find that the system is highly hysteretic: up to three different magnetic 
states can be simultaneously stable in one region of the parameter space (magnetic field and bias 
voltage) .We also discuss how the magneto-resistance can be used to provide additional information 
on the non-equilibrium peaks present in the nanomagnet spectroscopy experiments. 
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The observation of the current induced magnetic 
reversal^ in Ferromagnetic-Normal metal-Ferromagnetic 
(FNF) spin valves has been a real breakthrough in the 
field of spintronic. Indeed the spin torque effect at the 
root of this magnetic reversal allows one to control the 
magnetic configuration using electronic current only, i.e. 
avoiding the necessity of using magnetic fields. The po- 
tential for practical applications of this effect has been 
recognized very early and includes current driven Mag- 
netic Random Access Memories and Radio Frequency 
components. 

The theoretical description of this physics is usually 
done in two (relatively independent) steps where one 
treats separately the magnetic and transport degrees of 
freedom. One first considers the transport degrees of free- 
dom (electrons around the Fermi surface) and calculates 
for a given configuration of the magnetization the cur- 
rent and spin current flowing through the system. In a 
second step, one considers the magnetization dynamics. 
The divergence of the spin current is then identified as 
a source torque term for the dynamics of the magnetic 
degrees of freedom and is added to the equation of mo- 
tion of the magnetization. This physical picture, some- 
times associated with the sd model (the light s electrons 
are responsible for transport while the more massive d 
electrons carry the magnetization) is justified by the dif- 
ference of time scales associated with the fast transport 
degrees of freedom and the slow magnetic degrees of free- 
dom. In its simplest form, the magnetization dynamics 
of a mono domain ferromagnetic layer is described by 
the (phenomenological) Landau-Lifshitz-Gilbert (LLG) 
equation to which one adds the terms calculated from 
the study of the transport degrees of freedom 2 . 



The first corrections to this picture were considered re- 
cently 3 and lead to an enhanced magnetic damping con- 
stant due to the transport electrons, a feature that has 
been known experimentally for quite a long time-: the 
presence of non magnetic layers sandwiching a magnetic 
one induces a large increase of the damping measured in 
ferromagnetic magnetic resonance spectrum as spin cur- 
rent flows in and out of the magnetic layer. The aim of 
this paper is to study theoretically the simplest system 
where both magnetic and transport degrees of freedom 
can be studied on the same footing. Starting from a 
microscopic model containing only the original electrons 
we derive an effective equation for the magnetization dy- 
namics. The standard concepts of spintronics such as 
spin accumulation or spin torque appear within this for- 
malism and can be viewed as different aspects of the same 
physics. 

The simple system mentioned above is not a mere toy 
model of academic interest but is actually realized exper- 
imentally in tunneling experiments with ultrasmall mag- 
netic grains. Those nanomagnets (typically 3 — 4 nm 
thick grains of magnetic materials with a total spin 
So ~ 1000) are actual model systems for the study of 
the interplay between transport and magnetism. They 
are small enough so that their magnetization can be con- 
sidered as a large macrospin without spatial variations-. 
Indeed the magnetic exchange length is of the order of 7 
nm in cobalt so that degenerate magnons with non-zero 
momentum can be ignored. Moreover, their mean level 
spacing is of the order of lmeV, bigger than the tempera- 
ture so that the discreteness of the grain energy spectrum 
can be resolved^. A first set of experiments concentrates 
on the transport properties using single electron tunnel- 
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ing spectroscopy (connecting the grain to electrodes with 
tunnel junctions) to probe the grain&i. The differential 
conductance-voltage characteristic of such a device dis- 
plays peaks corresponding to excitations of the grain from 
which the grain's spectrum can be measured^. Another 
set of experiments focus on the magnetic degrees of free- 
dom of the grain by a direct measure of the magnetization 
using a microSQUID technique Reversal of magnetic 
moments as small as 1000 /is have been observed provid- 
ing information on the static magnetic properties of the 
system (the measure of the switching field as a function 
of the external field direction allows the reconstruction of 
the anisotropy tensor) and also on the dynamical mag- 
netic properties (switching times). 

The actual device considered in this paper is shown 
in Fig. and consists of a nanomagnet connected via 
tunneling barriers to magnetic contacts whose magne- 
tization can be non-collincar with the grain's. We use 
this system to illustrate various aspects of spintronics 
devices and in particular, current induced magnetic re- 
versal^, current induced spin waves excitationAi, spin ac- 
cumulation 12 , Tunneling magneto-resistancei&iii^ and 
magnetic damping through the conducting electrons^. A 
3 nm nanograin is not equivalent to a thin layer how- 
ever, and its smallness makes its physics more com- 
plex. We find that the presence of the Coulomb blockade 
phenomena™ strongly affects the spin dynamics. Its in- 
terplay with the magnetic degree of freedom leads to a 
"magnetic blockade" where electrons with a certain spin 
cannot enter the grain. We show that the phenomeno- 
logical LLG equation does not apply to the nanomagnet 
and has to be replaced by a new one ( Ea. l|24|l which to- 
gether with its corresponding phase diagram Fig.[S]is the 
main result of this paper) which incorporates the complex 
interplay between charge and spin dynamics. This equa- 
tion, that we refer as the NMD equation (Nano Magnet 
Dynamics), is not deduced from phenomenological con- 
siderations but rather derived from a microscopic quan- 
tum model. The derivation is interesting in itself as it 
sheds light on the origin of the time scale separation be- 
tween magnetic and transport degrees of freedom. Early 
results on the subject were obtained in Refiii, and some 
of the results presented in this paper were announced in 
a short publication^. 

The paper is organized in such a way as to provide two 
ways of reading. The reader not interested in technical 
details can jump directly to the introduction of section 
If VI where the NMD equation is presented and then to sec- 
tion ^where we study comparatively the phase diagrams 
of the LLG and NMD equation. In section we intro- 
duce our microscopic model for the nanomagnet. The 
magnet is described by a simple Stoncr like model and is 
connected to (non-collinear magnetic) leads via tunnel- 
ing barriers. An additional phenomenological coupling to 
an (bosonic) environment is introduced to correctly de- 
scribe the spin relaxation of the grain (via the coupling 
to phonons, spin-orbit coupling,...). In section^we de- 
rive the Master equation for the grain's charge and spin 
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FIG. 1: Left: schematic of the system. Several magnetic 
reservoirs of chemical potential fit and magnetization direc- 
tions rhi (in arbitrary directions) are connected to the mag- 
netic grain with tunneling barriers. The grain of magnetiza- 
tion m and easy axis z is also capacitively coupled to a gate at 
a potential V g . Right: The magnetization rh can be written 
in spherical coordinates (6, (f>) 



dynamics in close analogy to the one found in the stan- 
dard theory of Coulomb blockaded. In section, IIIII we 
apply this master equation to the Tunnelling Magneto- 
Resistance (TMR) and discuss how it could be used to 
get extra information on the spectroscopy of the grain. 
While the rest of this paper is devoted to magnetization 
dynamics, this section on charge dynamics is relatively 
self-contained and can be read independently by read- 
ers familiar with tunnelling spectroscopy of ultrasmall 
grains. Section IIVI is devoted to the derivation of the 
NMD equation in the semi-classical limit of a large grain's 
macrospin. We show that the NMD equation describes in 
a unified way the current induced torque and the mech- 
anism of magnetic relaxation through the metallic leads. 
Section [V] describes in details the phase diagram of the 
NMD equation together with the LLG equation for com- 
parison. 



I. MICROSCOPIC MODEL 

Our nanomagnet is weakly connected through tun- 
neling barriers to several magnetic electronic reservoirs, 
whose chemical potential, temperature and magnetiza- 
tion are respectively fii, Ti and mj, and is also capaci- 
tively coupled to a gate as shown in FigQ] Such a geom- 
etry corresponds to the experimental geometry of Ref. 7. 
We neglect any spatial variation of the magnetization 
in the grain due to the large magnetic exchange length 
(compared to the size of the nanomagnet). In order to 
describe properly the relaxation properties of the grain 
(i.e. those not due to electronic transport through the 
leads), we introduce a second coupling term to an en- 
vironment described by bosonic degrees of freedom (not 
displayed in Fig^}. This bosonic bath represents all the 
intrinsic degrees of freedom that can provide spin re- 
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laxation (phonons, non uniform magnons, impurities...). 
As we shall see the presence of this term will lead to 
a Gilbert type damping term in the magnetization dy- 
namics. Other types of relaxation mechanism that con- 
serve spin have not been included although their pres- 
ence might modify our results quantitatively (they affect 
the charge dynamics). The presence of non-equilibrium 
peaks in the spectroscopy experiments^ is however a 
strong indication that charge relaxation rate is fairly low 
in those systems. 



A. Isolated magnetic grain 

We start our discussion by introducing a simple Hamil- 
tonian that describe the isolated nanomagnet with an 
easy axis along the z axis. It is similar to the one dis- 
cussed ir&iS, and consists of a charge (c) and spin (s) 
part, 

H grain = ^ grain "f" ^graim (l a ) 
#g C rai„ = E ^S aa d aa + E C (N - N g )\ (lb) 



-^grain — JS ■ S — - S ' z Tv~fHS z 
DO 



(lc) 



Here, d' aa (d ar7 ) creates (destroys) an electron on a level 
a with energy e a and spin a along the z-axis. N is the 
number of electrons in the grain, Ec the charging energy 
and N g is proportional to the gate voltage. The exchange 
energy (oc J), uniaxial anisotropy (oc ft) and Zeeman cou- 
pling (oc H, we consider only a magnetic field along the 
z axis) are function of the total spin S, whose expression 
in term of the electrons present in the system reads : 



dacri&<ri<T2 

a,tTi ,(T2 



dac 



(2) 



For k — and J smaller than the mean level spacing at 
the Fermi level, this model is known as the "Universal 
Hamiltonian"^SiSi and has been widely used to describe 
quantum dots on energy scales below the Thouless en- 
ergy. It corresponds to the limit of zero dimensions (i.e. 
where the physical quantities are dimension independent) 
which is obtained when all the time scales involved are 
larger than the time to explore the entire volume of the 
grain. When J is increased, the system undergoes a fer- 
romagnetic (Stoner) instability and acquires a finite spin 
So- The eigenstates \A) of iJ gra in of energy Ea are read- 
ily found by observing that the charge and spin parts of 
the Hamiltonian commute with each other so that the 
number n a = d) aa d a(J ) = 0, 1, 2 of electrons in level 
a, the total spin S, and the spin S z along the easy axis 
z of the dot are good quantum numbers. 

The analysis of the low energy spectrum of i/ gra in 
has been done extensively in&iS. Here, we reproduce 
the main results for the sake of clarity, (i) The charge 
part of the ground state is shown in Fig[2j (a): N s = 



J2 a n a(2 — n a ) one-body levels are singly occupied in or- 
der for the system to have a non zero spin So = N s /2. 
The Fermi level is then split into majority and minority 
electrons, with different mean level spacings, respectively 
Am and A m . (ii) It costs a huge amount of energy cx J So 
to switch the spin of one electron without changing the 
filling factors (n a ) of the one-body levels. J So is of the 
order of the Curie temperature (around 1000K) which is 
very large compare to all other energy scales considered 
here. Hence we can safely ignore all the excited states 
with S < Ng/2. We thus assume, for each eigenstate 
that S — N s /2. For these low energy excited states, 
there is no additional degeneracy, and an eigenstate is 
entirely described by the set of n a and S z : 



\A) = \{n a },S z ) 



(3) 



(iii) there are three types of low energy excitations, as 
sketched in Fig|2 (a) and (b) . Following-, we shall call 



(a) 




(b) 



Majority minority 



FIG. 2: (a): Occupation of the different one-body level in the 
ground state of the system. Squares (resp. circles) represent 
the majority (minority) electrons with mean level spacing of 
Am (A m ). A,B,C and D represent low energy excitations, see 
text, (b): idem but the exchange energy has been taken into 
account in addition to the one body energy to show explicitly 
that the spin accumulation excitation (C and D) have energies 
of the same order as the particle-hole excitations (A and B). 

them particle-hole excitations for majority (A) and mi- 
nority (B) electrons, spin accumulation (C and D: a par- 
ticle jump from a majority (minority) level to a minority 
(majority) level hence changing S by one unit) and spin 
precession excitation (E, not shown in the figure: n a and 
S are unchanged but S z is changed). We note that the 
huge amount of one-body energy (cx J So) that the D ex- 
citation costs is almost exactly canceled by the gain in 
exchange energy, so that finally all these low energy ex- 
citations are of the same order J ~ Am ~ A m <C J So. 
The spin precession excitation have an energy of order 
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k. (iv) the mean level spacing and exchange energy scale 
as the inverse of the grains' volume while the anisotropy 
parameter k is volume independent. For the grains con- 
sidered in Refs 6, 7. (3nm, around 1000 atoms) however, 
both have the same order of magnitude (k w 2QQmK 
and A m < Am of a few Kelvin). 

We note that although we start with only electronic 
degrees of freedom, we end up with low excited states 
of charge type (A,B,C, and D) and magnetic type (E). 
The dynamical interplay between the two will be induced 
by the relaxation channel (coupling to the bosonic bath) 
and tunneling coupling to the leads. 

B. Coupling the nanomagnet to the leads and the 
environment 

Once the isolated grain is understood, we introduce 
the coupling to the leads and to the environment. The 
full Hamiltonian of the systems reads : 

H = //grain + #leads + H t + H cnv + H g (4a) 

H t = ^2 l ak a ^ ak(J d aa + h.c. (4b) 

akaa 

H\c&d = ^2 Zk@V ak pCakp (4c) 
ak/3 

H s = g B (S+ < t> + S-^) (4d) 

In this expression, c aka (c a ka) create (destroys) an elec- 
tron in lead a with momentum k and spin a along the 
z-axis. The (possibly) ferromagnetic leads are described 
at the mean held level, using different energies ek/3 for mi- 
nority and majority electrons. Here, the index f3 = M, m 
(Majorityminority) is understood along the actual mag- 
netization of the lead, which is not necessarily collinear 
with the easy axis of the grain (z). The c^ fc(T are re- 
lated to the c' ak g (that creates electrons along the mag- 
netization of the leads) through a SU{2) rotation matrix 
c lka = Y^nKAkfi with Ra = exp[-*0«c? • n a /2}. 9 a is 



the angle between the magnetization m a and the z axis, 
while ri a = zx m a /\zx m a \. We assume that no spin-flip 
takes place during the tunneling processes. 

H env describes the environment in term of a bosonic 
field (f>. It is an effective, phenomenological description 
of various intrinsic relaxation processes of the dot, like 
spin-orbit scattering or phonon-magnon scattering that 
couple the grain's spin to bosonic degrees of freedom (like 
phonons or non-uniform magnons). A detailed study of 
these processes is beyond the scope of this paper. We 
restrict ourself to a phenomenological description of in- 
trinsic damping mechanisms through spectral density of 
the boson </>, which will be denoted by Pb(w)- We assume 
that pB is wide and slowly varying on the frequencies 
scales at which the spin precesses so that the environ- 
ment do not introduce new time scales in the problem. 
Moreover, we restrict ourselves to a bath coupled to the 
spin part only, with a coupling constant g B . We do not 
introduce a coupling to S z (which do not provide relax- 
ation and hence would play no role in the following) or 
to the charge degree of freedom (we suppose that the 
corresponding relaxation time is long as indicated by the 
presence of non-equilibrium peaks in the tunneling ex- 
perimentsi). 



II. MASTER EQUATION. 

The coupling to the leads and environment induces the 
dynamics of the grain's degrees of freedom. Here, we use 
a direct extension^ of the theory of Coulomb blockade in 
the sequential tunnelingi^ regime. The dynamics is de- 
scribed by the Master equation for the probability P A (t) 
for the system to be in state \A) at time t. This equa- 
tion has a simple physical interpretation: the different 
(tunneling and interacting with the environment) events 
occur incoherently (with rates oc T° and oc respec- 
tively) and induce a random walk for the state of the 
system. The master equation is the equation of conser- 
vation of probability for this process. It reads, 



dP A 

dt 



E j r aa|<^l<U^>| 2 \n a (AE)P A , (l-n a (AE))p A 



+ 1^,1(^14,1^)1 



-n a (-AE)P A 



(l-n (-A£?))fV |+7r<?i E e\(A'\S x -ieS y \A)f 



A', e=± 



PB {eAE)n B (AE)P A , + p B (eAE)n B (-AE)P A 



(5) 



where AE — E A ~ E A i is the energy difference be- 
tween state \A) and \A'), n a {E) = 1/(1 + e( E "^)/ feT ) 
(n B (E) = l/(e E / kT - 1)) is the Fermi (Bose) function, 
T the temperature and p a the chemical potential of lead 



a. r^ Q is the tunneling rate (see Eq. (|T4"}) and Eq. lfTS)l 

for the precise definition) . The term proportional to T° a 
in Eq.(jSJ corresponds to the usual master equation for 
tunneling processes, while the second term oc T B is the 
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corresponding term for bosonic degrees of freedom. 

As we shall see, the validity of the master equation 
approach imposes that the time between two different 
events {H/T^ a and H/Tb) is large compared to the coher- 
ence time of the leads/environment h/(kT). The exten- 
sion to a magnetic grain is not straightforward however 
as the presence of non-collinear magnetic contacts could 
lead to off-diagonal terms in the density matrix (there is 
no good quantization axis for the spin) . We find that the 
master equation remains valid when internal precession 
is fast enough to average out contributions transverse to 
the easy axis. The time needed for S z to change signif- 
icantly is of order Sq/T (one tunneling event can only 
change S z by half a unit) so that this condition implies 
r -C kSo/H. 

In this section, we perform two tasks. First, we give 
a proper derivation of Eq.JSJ) and discuss its domain of 
validity. Secondly, we evaluate it for the Hamiltonian of 
our system Eq. l|4a|l and arrive at the specific form for the 
magnetic grain, equations Q), i|TT?|) and QH). 



A. Quantum master equation 

Our starting point for the dynamics of the system is the 
quantum master equation for the reduced density matrix 
p of the dot, as derived for instance in2£ : a generic sys- 
tem a (the magnetic grain) coupled to an environment 
b (the leads and the bosonic environment) is described 



by a Hamiltonian Htot = H a + gijAiBj + Hj, where the 
operators A and B act on a and b respectively. The 
quantum master equation for the reduced density matrix 
p (where the environment degrees of freedom have been 
traced out) is obtained by writing the evolution equation 
of the full density matrix in second order in the couplings 
gij . The trace on the environment can then be taken and 
one get an effective equation for the system evolution, 

dtp = -C9ij9ki [ dr [Mt),A k (T)p(r)] (Bj(t)Bi(r)) 



+ [p(r)Mr),A k (t)}(B j (r)B l (t)) (6) 

where C = 1, — 1 if the operators A and B commute or 
anticommute. If the environment is much larger than the 
system, we can assume that it is in thermal equilibrium 
and calculate the correlation functions for the Bj opera- 
tors accordingly (note that this is not true for the system 
since it is microscopic and its out of equilibrium steady 
state density matrix is determined by its couplings to the 
environment, even at order zero in those couplings 23 ). In 
the following, (• • •) stands for the thermal average of the 
lead. Given the form of the microscopic Hamiltonian Eq. 
(|laj) . at second order in the couplings, the master equa- 
tion has the form : 

dtp = £>Leads/0 + T> cnv p (7) 
with X>Loads/6 ~ 0(t 2 ) and T> cnv p ~ 0(g 2 B ) and we can 
compute the effect of the electronic leads and the bosonic 
bath separately. 



B. Coupling to the leads 

First, let us derive the first term T^eads/^ of Q). Ex- 
panding Eq. ijBJ in the tunneling matrix elements, we 
obtain : 



dr 



d a At),dl al (T)p(r)} A«'/(t-r) 



olol 1 aa 

+ \p{r)dlAr),d aa {t) B a J\t-T) + h.c.\ (8) 



with the definitions, 



Ka'it-r) = l^t akoi t* aka ' (cak,'(0)cl k<J (t-T 



(9a) 



Caka{t) 



(9b) 

(9c) 
(9d) 



Eq. in general is rather complicated since the deriva- 
tive of p at time t involves the density matrix at former 
times. Hence the system has "memory" . The leads corre- 
lation functions (t) and (t) can be calculated 
and behaves typically like 



Ala it) OC 



constant t < h/ep 
l/t h/e F < t < k B T 



-■Kk B Tt/h 



(10) 



* > k B T 



where T is the temperature of the reservoirs and cf the 
Fermi energy. Therefore this "memory" has a correla- 
tion time of h/ksT. If the density evolves on time scales 
longer than that, p(r) can safely be replaced by p(t) in 
Eq. (|SJ • This is known as the Markovian approximation 22 
and is valid for H/kT <C H/T where T is the finite 
width of the grains' levels due to the coupling to the 
leads(calculated using the Fermi Golden Rule, see the for- 
mal definition in Eq. 1151) . The domain of validity of the 
Markovian approximation (fcsT 3> T) is consistent with 
the fact that Eq. (JHJ is obtained by formally expanding in 
the hopping amplitude, and hence T should remain the 
smallest energy scale of the problem. In the other limit, 
k B T <C r, we have A^J (t) oc l/t, the integration over 
r leads to a divergence and Eq.© breaks down entirely. 
Once the Markovian approximation has been done, we 
project Eq.JSJ on the many-body eigenstates of the iso- 
lated dot, and carry out the integration over t. Noting 
Pab = (A\ p \B) the projection of the density matrix, we 
obtain 
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T^Lc&dsPAB = ~ ^ EA ~ ED)t pDB (A\ d aa \C) (C\ d^, \D) (E c - E D ) 

- ^ Ec ' EB)t PAc (C\ <C \D) (D\ d ao \B) BZ{E C - E D ) 

+ e ^ A -E c+ E D -E B )t f)cD {A] | C) (jD | ^ ]B) AZ \ Ea _ Ec ) 

+ e ^ A -E c+ E D -E B )t~ pcD (A | ^ ^ (D | ^ {Ed _ Eb) + {a ^ B y (n) 



where we introduced 

poo 

A(E) = dt e- lEt A{t) (12) 
Jo 

Ea. (|ll|l reduces to the (first part of the) master equa- 
tion Eq. JSJ) if we ignore completely the off diagonal part 
of the density matrix and note Pa = SabPab its diag- 
onal part. Such an approximation, although very com- 
mon in Coulomb blockade theory needs to be justified 
in the present context. Indeed, when the magnetizations 
of the reservoirs are non-collinear with the easy axis of 
the grain, neglecting the off diagonal part of the density 
matrix amongst to neglect all the dynamics in the xy 
plane (as the matrix elements paa depend only on the z 
component of the magnetization). In other words, in this 
approximation an electron tunneling along the x axis is 
equivalent to an electron tunneling along the y axis: if 
two reservoirs have their magnetization lying in the xy 
plane, then the physics is independent of the relative an- 
gle between their two magnetizations. A case where the 
off-diagonal elements of the density matrix have to be 
considered has been discussed irj24. 

On a qualitative level, we find that the off-diagonal 
part of the density matrix oscillates at a frequency set by 
the excitation energies (the phases in Ea . (|1 ip ^ while the 



density matrix evolve on a rate T. Hence, when the exci- 
tation energies are larger than T, the off-diagonal density 
matrix averages to zero. For our model, this translates 
into k>1. In this limit, which is achieved for instance in 
the grains studied in^i, the magnetization of the grain 
precesses several times around its easy axis (z) in be- 
tween two tunneling events, so that the transverse de- 
grees of freedom (along xy) are averaged out, and the 
only relevant quantity is the projection of the magneti- 
zation along the z axis. In this limit, all the physics of 
non-collinear reservoirs can be understood as if the mag- 
netizations of the reservoirs were all lying on the z-axis. 
This argument can be made more precise as explained in 
Appendix A. The master equation is of course also valid 
when the leads magnetization are parallel to the easy axis 
z. 

The above argument is rather general, and does not 
depend on the particular form of -f/g ra in- As we shall 
see below, the magnetic grain displays a separation be- 
tween the charge and spin relaxation times (respectively 
of order 1/T and So/T). This shall increase the range of 
validity of our approach to k> T/So (see below section 
IE). 

Eq.JITJ reduces to the first part of Eq.0, 



C.aacr 



V Lcads p A = \{A\d atJ \C)\ 2 (pcna{E c - E A ) - p A n a {E c - E A )) 

+ | (A\ 4„ \C) \ 2 (p c n a {E A - E c ) - p A n a (E A - E c j) 



(13) 



where n a {x) = 1— n a (x) and T° a is the effective tunneling 
rate from reservoir a on the state a along the z-axis. It is 
connected to the tunneling rate along the magnetization 
direction rfi a of reservoir a through, 

rL = cos 2 ^rL + sin 2 |rf Q 

rL = cos 2 ^rf Q + sin 2 ^rl (w) 

where 8 a is the angle between the z axis and fn a , while 



^aa {P — U, D) is given by the Fermi Golden Rule: 

k 

Eq. (|13|) was used previously for reservoirs whose magne- 
tization was parallel to the easy axis of the grain. Here its 
regime of validity is extended to non-collinear reservoirs. 

The average electrical current flowing from reservoir a 
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onto the grain reads, 



A,C,ota 



\(A\di a \C)\ 2 n a (E A ~E c ) 



spin part is a Clebsch-Gordon coefficient and reads (Cf. 
Appendix C, formulas C.25 and C.26 o££), 



5 , S z 



dt 



5, 5, 



-|(A|d aCT |C)| 2 n Q (£;c-^) 



(16) 



We now turn to evaluate Eg. If 3fl specifically for the 
case considered in this paper. We restrict ourself to the 
Coulomb blockade regime where Ec is large enough so 
that only states with N and N + 1 particles need to be 
considered, and we note V g the corresponding charging 
energy (which is tunable with the gate voltage). Let us 
define p + ({n'p}, S' z ) (resp. p_({np}, S z ) the density ma- 
trix in the sector with N + 1 (resp. N) particles with 
S' = \T,pn' p {2-n' p ) (resp. S = \ £^(2-^)). The 

matrix elements, (S z , {np}\ S aa S' z ,{n'p}^j factorize into 

a charge part and a spin part. The charge part, noted 
Fa{{np\,{n'p}) is equal to zero except when np = n'p for 
every (3 ^ a and n' a — n a — 1 where its value is one. The we finally obtain 



z - — —8s>,s+3?-8s> s ,s,+% (17) 



where rj a = +f (rj a = — 1) for a majority (minority) 
state. After calculating the energy differences (the one 
body energies are counted from their respective minor- 
ity/majority Fermi level), 

AE = e a + V g — J^±^ 



AE' = e a + V q - J- 



K 


S z a - 


K 


"*5 


Sq 












(18a) 


K 


S' z a + 


K 










So 




45^ 


2 








(18b) 



V^ s p + ({n' },S' z ) = Yl K a F a ({n p },W }) S '^' hr ' S '^' 



25' + 1 - j] a 



v a {AE')p. (W, ^ - |) - na{AE')p + (W p }, 



©Lead 8 P_({^},^)= £ KaFa({n0},{n'p})- + ''" fTS 



aaa{n' 3 } 



25 + 1 



i a {AE)p + ({n^}, 5 2 + |) - n a (AE)p- ({n p } : S z ) 



(19a) 



(19b) 



It is important to realize that four types of tunneling 
events can occur: The electrons can tunnel either on a 
majority or minority state, with a spin either up or down 
(with respect to the easy axis z). This apparently para- 
doxal statement just reflects the fact that we measure 
the spin along the easy axis and not along the magneti- 
zation of the grain. When the system is near equilibrium, 
the Clebsch-Gordan coefficients almost vanish when a up 
(down) electron tunnel on a minority (majority) state, 
making the distinction between minority/majority and 
up/down electrons meaningless. However, when the sys- 
tem is put out of equilibrium (as it will be in the presence 
of a bias voltage) such a distinction is relevant. 

C. Coupling to the bosonic environment 

In this paragraph, we derive the second term T> cnv p of 
(|7|). The quantum master equation for the environment 



reads using Eq.© and noting = S x ± S y : 



d t p=-9% / dr S+(t),S-(r)p(r) (#0)0* (r-t)) 

J — OO 

-g% f dr\s-{t),S + {T)p{r)](^{<S)ct>{r-t)) + h.c. 

J — OO 

(20) 



As in the previous calculation, we use the Markov ap- 
proximation, neglect the off diagonal terms of the density 
matrix projected on the eigenstates of the grain. We ob- 
tain the second part of Eq|5] The specific form for our 
problem, reads 
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■D cnv p ± ({n a }, S z ) = g% l\(S z \ S+ \S Z - 1)| \A^( E s.-i ~ E s Jp ± ({n a }, S z - 1) - B^E s „-i ~ E s Jp±({n a }, S z ) 
+ \(S Z + 1| S + \S Z )\ 2 (b^Es; - E Sz+ i)p ± ({n a }, S z + 1)- A^(E Sz - E Sz+ i)p±{{n a }, S z ) 



r 



where we used the definitions 



ME) 
B*(E) 



-hS 



-S 



So 



<^(u)^(0)) 



,-iEt, 



(22a) 
(22b) 
(22c) 



The imaginary part of A$ and only renormalizes the 
precession frequencies of the spin— i, so we will drop them. 
After standard manipulations, we reduce (|22a|l to 



ReA^{E) = wn B (-E)p B (-E) (23a) 
ReB^E) = -Trn B (-E)p B (-E) (23b) 

where n B (E) = l/(e E ^ kT — 1) is the Bose function and 
h B (E) = n B (—E). Note that as the bosonic environment 
couples only to the spin degree of freedom, it is diagonal 
in p + , p_ . Equations I|7I19I21|I define the master equation 
that will be used in the rest of this paper. 



III. GRAIN'S SPECTROSCOPY AND 
TUNNELLING MAGNETO RESISTANCE 



One important aspect of the Master equation (|7I19I21|) 

is its validity for non-collinear systems. Indeed in real 
samples, it is very difficult to have the easy axis of 
the grain aligned with the magnetization of the contact. 
However, Ea. (|14|) can be interpreted as if the anisotropy 
axis of the grain were aligned with the magnetization di- 
rection of the contact. Hence, the tunnelling magneto- 
resistance can be measured (in a sample with one mag- 
netic lead) by switching the direction of the lead magne- 
tization while keeping the magnetization of the grain un- 
changed (this could be done for instance using permalloy, 
whose coercitivity of the order of 10 mT is much smaller 
than the anisotropy of the grains, typically lOOmT) irre- 
spectively of the angle between the two magnetizations. 
As we shall see, those TMR experiments can provide 
extra information on the spectroscopy of the grain and 
in particular should help to differentiate non-equilibrium 
peaks associated to spin accumulation excitation from 
those associated to particle-hole excitations. 

The current- voltage characteristics of a magnetic grain 
with non magnetic lead has been analyzed extensively 



both theoretically 8 * 9 - and experimentally 6 ^. The differ- 
ential conductance dl/dVbias shows peaks that allow one 
to extract information about the low energy spectrum of 
the grain (assuming without loss of generality that upon 
applying a bias voltage Vbias across the sample the chem- 
ical potential of lead 2 is unchanged while the one of lead 
1 is increased by eVbi as )- The basic rules to get a peak at 
a V£ ias are (i) there must be two eigenstates of the grain 
|Ajv) and |Ajv + i} with respectively N and N+l particles, 



such that Ea, 



E An = el/ b * ias (being in \A N ) the sys- 



tem must be able to get enough energy from the lead to 
reach | j4jv+i ) in a tunneling event), (ii) There must be a 
non zero matrix element (AN + i\d) aa \Ais[) . (iii) The prob- 
ability of being in state \Am) must be non zero. The case 
where \Am) is the ground state of the grain with N elec- 
trons is the simplest to analyze and has been used to ex- 
tract the one-body energies e Q of non magnetic grains^ 7 .. 
All the other peaks are referred to as non-equilibrium 
excitations 2829 . In non magnetic grains, particle- hole 
excitations are at the origin of extra peaks, provided 
the charging energy has mesoscopic fluctuations, i.e. de- 
pends slightly on |^4jv). In addition it was found that in 
magnetic grains, spin accumulation excitations also give 
peaks while spin precession excitations do not. 

In order to illustrate the interest of the TMR experi- 
ments, we have calculated numerically the exact solution 
of Eq. (|7I19|) without bath and plot in FigOl the differen- 
tial conductance (left panel) and the difference between 
the current in the aligned and anti-aligned configuration 
(right panel) as a function of the bias voltage Vbi as and 
gate voltages V g for an illustrative sample. The station- 
ary solution of Eq. J7J amounts to solve a linear equation 
of the form AX = which is efficiently done numerically 
by separating the diagonal D and off-diagonal part E 
of A = D — E and multiplying iteratively a initial vec- 
tor X by the matrix D~ X E until convergence (which is 
guaranteed by the fact that A is a stochastic matrbi^i). 
In the left panel for V g < a w 2 we see a typical spectra 
where each peak (a,c,and e) is associated to a different (in 
this case minority) level (cartoon 1). For V s > a, (car- 
toon 2 and 3) we see that two extra "non-equilibrium" 
peaks appear (b and d). These peaks are associated to 
spin-accumulation excitations and have been described 
in detail inS (an electron enters on the minority level (a) 
leaves from a majority level which decreases the exchange 
energy of the grain hence the bias voltage needed for an 
electron to tunnel on the (c) minority level) . On the right 
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panel, we see a similar structure (except that total cur- 
rents have been plotted not conductances) with an addi- 
tional structure at V g < (3 « 3 which is barely visible on 
the left panel. This additional structure comes from the 
fact that there can be more than one spin-accumulation 
excitation that gives a peak approximately at the same 
bias voltage (in our example the electrons can leave from 
two different majority levels, cartoon 4) and while the 
positions of the two peaks might be hard to resolve, their 
contribution to the magneto-resistance can be very dif- 
ferent. We also note that the sign of the TMR shown in 
FigElis negative almost everywhere (the tunneling hap- 
pens on minority levels) except for the upper right corner 
V g > (3 and Vbi as > e where it becomes positive. Hence 
the sign of the TMR can be controlled, at least for the 
choices of parameters used in this example, by external 
voltages. 
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FIG. 3: (Color Online). Differential conductance dl/dVbins 
(left panel) and difference between the current in the aligned 
and anti-aligned configuration 1(0 = 0) — 1(0 = n) (right 
panel) as a function of the bias voltage Vbias and gate voltages 
Vg for an illustrative sample. The magnitudes (in arbitrary 
unit) increase when going from blue to red. The cartoons 
represent the different type of tunneling events occurring for 
different values of Vg and Vbi as - The levels on the left of each 
cartoon represent the minority levels while the majority are 
on the right. The dotted line stands for the Fermi level of 
the right (non magnetic) contact while the electrons entering 
from the left (magnetic) contact have an energy which is Vbias 
higher. 



IV. MACROSCOPIC SPIN LIMIT 

In actual nanomagnets, the total spin is large (around 
1000) which considerably simplifies the master equation 
into a differential equation for the magnetization dynam- 
ics analogous to the LLG equation. In this section, we 
perform this derivation. First, we take the large spin 
Sq limit in the master equation, and replace the discrete 



degree of freedom S z (projection of the spin along the 
easy axis z) by its continuum version u ~ S z /Sq. Sec- 
ond, we find that charge and spin equilibrium occurs on 
two different time rates, (fast T for the charge and slow 
T/So for the spin) so that we can integrate the charge 
dynamics. The result has the form of a Fokker-Planck 
equation for u. The associated Langevin equation (the 
NMD equation) is the chief result of this paper and reads, 



R(u) 



l2D(u) 



R(u) = (1 



as(h + 2ku) 



(l-q)h(u)-(l + q)f l (u) 



D(u) = (1 - u z ) 



AS Q {l + qriu) 

h + 2ku 



(24a) 



(24b) 



a B 



ZU1 l 2kT I 



(l-gu 2 )/ T (u) + (l + gM 2 )/ i ( M ) 

8S (l + qvu) 
u2(/ T ( u )-/ i ( u )) 2 -(/ T ( u ) + / i ( u )) 2 ^ 



8S T(l + qr]u) 



(24c) 



where we have introduced the following notations: the 
magnetic field along z is h = "/H, £ t is a white noise 
(£t£f) = ~ an d the Gilbert damping term olb 
is the intrinsic damping of the isolated grain. r\ = ±1 
depending on whither the tunneling occurs on a major- 
ity/minority state. The functions /o-(m) (with a = ±1 for 
up/down spin) contain all the dependence in temperature 
T and leads chemical potentials /i a , 



f<r{u) =^r>(e-/tau-/i| 



AM 



(25) 



where the offset energy e includes the charging energy 
for an electron to enter the grain, n is the uniaxial 
anisotropy and n(E) is the Fermi function. For simplic- 
ity, we have restricted ourselves to the case where the 
tunnelling events occur on only one one-body level with 
tunneling rate T^. q is the spin-torque asymmetry (or 
the average polarization of the tunneling rates) and T 
the total tunneling rate, 



^2a rl + ri 



ri 



(26) 



At the classical level, the z component of the spin has a 
non-linear deterministic evolution, corrected by a small 
stochastic noise. The deterministic part takes a form 
similar to the LLG equation (See below Ea . Q44JI 1 ■ except 
from the current induced term which is modified by the 
Coulomb blockade phenomena. The analysis of Ea. l|24() 
will be done in section Ivl while the rest of this section is 
devoted its derivation. Hereafter, Equation (|24|l will be 
referred to as NMD (Nano Magnet Dynamics). 
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A. Large spin limit 

As we restrict ourselves to tunneling events on one 
orbital a, we drop the index a. The system can be 
in one of two spin states with spin Sq (N particles) 
or 5q = So + r)/2 (N + 1 particles). For later conve- 
nience we perform the large spin expansion in 1/ S m with 
S m = Sq + 77/4, and denote k/S m = k/Sq. Moreover, we 
assume g 2 B = 0(T/S m ) so that the intrinsic spin relax- 
ation rate will be comparable to the relaxation induced 
by the presence of the contacts. In the classical limit, 



the magnetizations are continuous variables — 1 < u < 1, 
such that S' z = S' u and S z — Squ. We denote by p± 
the functions (of u) such that (normalization factors are 
chosen for convenience): 

p + (S z ) = ^LpJ§) , p_(S.) = $Lp-(%) (27) 

Introducing the notation x a = e + V g — J 1 ~ l 2 2 '' — kug — 
ha/2, the master equation I|7I19I21[I now takes the form, 



d t p + (u)=^T a ^-^ 

a,cr 

+ V cnv p+(u) 



n a [x a + (1 - rjau) ) —p 



Sq V Sq 



(28a) 



25*o + l 



S Q /5 m+| 



45, 



■(1 - ■qau))p-{u) 



+ X>envP- («) 



(28b) 



We now proceed with the expansion in l/<S m which assumes that that the functions p varies slowly with u, such that 
\d u p±(u) \ < |p±(w)|5 m . We get, 



d t p + {u) =2J -y- (1 + W u ) ^a(^o-)p_(w) - n a (a; CT )p + (u)J + O {T/S m ) 

a,cr 

d tP (u) = - 9„( Jlcads + j env ) + o (r/s£) 



(29) 
(30) 



jcnv = - ng%pB(-h - 2ku) 



Jlcads — 



1 \ f 1 - M 2 ) 

1 - u 2 + —)p{u) + K ' (1 + 2n B (-h - 2Ku))d u p(u) 

a(l — r/au) 2 



4&, 



fn a (a; CT )p_(u) - h a (x a )p + (u) 



+ 7^— cr(l - M 2 )(n a (a;<T)jO+(M) + (2?y - l)n a (x o -) / 0_ (u)) 
— — (1 - rjau)(l - u 2 )(n a {x a )d u p+(u) + n a {x a )d u p-(u) 

4:0™ V 



(31) 



(32) 



where we have introduced p(u) = p+(u) + p-(u), the to u = ±1, see below) the charge degree of freedom is at 
reduced density matrix for the spin only, traced over the equilibrium and we have, (up to terms of order T/S^) 
charge degrees of freedom. 



B. Spin dynamics 

We note that the terms of order T only couple to the 
charge degrees of freedom. The terms that affects the 
spin (it) are of order T/S m so that in leading order, (and 
we will be interested in the sub-dominant order only close 



2 a(u) + r)ub{u) 
P+W = =; r— P(u) (33a) 

1 1 + qi]u 

2 a(u) + ryub{u) 

P-( u ) = r — rr p( u > ( 33b ) 

1 1 + qrju 
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where 

per per 

a(u) = -fn a (x a ) a(u) = J2^-h a (x a ) (34a) 

a. a a,er 

per per 

b(u) ee ^cr-^-n a (x CT ) 6(m) ee y^g-j-n a (x tT ) 

a.cr a,cr 

(34b) 

(we recall that n(x) ee n.(— x)). Using the charge equilib- 
rium relations, Ea. (|30|l finally takes the form of a Fokker- 
Planck equation and reads, 



d t p{u) 
R(u) 

D(u) 



= d v 



R{u)p{u) + -^-d u [ D(u)p(u) 



a 



a(u)b(u) — a(u)b(u) 



(35a) 



TS m x 1 + qr\u 

ng%(l - u 2 )p B {-h - 2ku) + O (T/Sl) (35b) 
1 



2TS 7 



,„ 9 . a(u)a(u) — u 2 b(u)b(u) 

-(1 - u) 

1 + qrju 



1 + 2n B (-h - 2ku)) + {T/Sl) 



(35c) 



where i? and D have be calculated only at dominant or- 
der. The sub dominant correction to R will be stud- 
ied below. Equation (|35a|l is in the Ito formal, and the 
corresponding stochastic differential equation can there- 
fore be identified straightforwardly. Denoting as = 
Trg 2 3 (dpB(0)) > (we assumed that is slowly vary- 
ing on the frequencies scales at which the spin precesses 
to expand ps linearly), this stochastic differential equa- 
tion is (at dominant order S m = So) our central result 
Eq.l|24p. Note that we have neglected corrections to R(u) 
of order 0(T/S^). 

The electrical current flowing through the system 
(from lead a to the dot) can also be calculated within 
the continuous limit from Eq. Ijltjfl : 



Ia = e J' du PM £ J2 (1 + (1 + W ' U) ^ 



l 



rjpu 



Tin x 



(n a (x<r) - n b (x' a fj (36) 



At a stable fixed point u, the density matrix is a Dirac 
peak and we get : 



j = l^^[ 1 +1 M )( 1 +W'«) r a rJ ' 



aa' b^a 



l 



rjpu 



(n(e — nan — h— — p, a ) — n(e — nau — h— — pb) 

(37) 

At this point let us rediscuss the range of validity of 
the NMD equation. In the derivation given in section ITT1 



of the master equation, we have used the condition that 
the tunneling rates are small compared to the excitation 
energies of the grain (to neglect the off-diagonal terms of 
the density matrix). For the grain, this condition trans- 
lates into it » T. The derivation was general, and did 
not make use of the particular structure of the grain's 
Hamiltonian. In particular, as we have seen, the charac- 
teristic time scale of the charge and spin dynamics differ 
by a factor 1/Sq. This extends the validity of the NMD 
equation to k > T/Sq- Physically this can be seen as 
follows: as the spin So is large and a tunneling electron 
can only change it by 1/2, the time needed for S z to 
change significantly is of order So/T while the time for 
the internal precession of the magnetization is of order 
1/k. Hence, when k 3> T/So, the system has made many 
precessions before the spin had the time to change signif- 
icantly and the dynamics of transverse magnetization is 
therefore averaged out. However, it would be necessary 
to take the transverse dynamics into account to compute 
the 1 /So corrections to the electrical current Eq. (|37|l . 

The noise term in Eq. (|24fl is of order 1/Sq smaller 
than the main deterministic term and do not play a ma- 
jor role in general. There are regimes however (when 
the bias voltage is high enough) where the deterministic 
term completely vanishes and the dynamics of the mag- 
netization becomes completely diffusive. Such a regime 
has been studied in detailed in Refiii. In order to make 
contact with this work, let us note that there are sub- 
dominant order corrections to R(u), which we have ne- 
glected so far. When u — > ±1 however, R(u) vanishes and 
the terms of order T/S^ become important. To make 
contact with the results o£i£, we derive the evolution of 
the spin close to the "pole" u = 1. The correction SR(u) 
of order T / Sf n to R(u) reads : 



5R(u) 



1 



TSl 



db — ab + u(bb — da) 
1 + prju 



+ (l-u 2 )f(u) 



where / is a regular function in u = ±1. For u close 
to 1, we linearize with u = 1 — A4/S and look for the 
(linear) equation of evolution of Ai. In order to ob- 
tain an equation for the (noise-) average (A4), we take 
advantage of the Ito form of the Fokker-Planck equation 
and the associated Langevin equation (|24|l . which implies 
{D(u)£ t } = 0. Putting together the contributions com- 
ing from Ea. (|24|l and Eq. l|38|l . we obtain (at the order 

r/^ 2 „), 



j t {M) = 1+ -{ 1 -- 1+ )(M) 

_ £arcq(sT)rl£ 6 "6fo)r£ 
r(i + p»7) 

Y, a na{x ] )TlY, b n b (x i )Ti 

+ " ra + pT?) 

which is precisely the equation obtained in 17 . 



(39a) 
(39b) 

(39c) 
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V. PHASE DIAGRAM 

We now turn to the study of the phase diagram of 
the NMD equation 124|l . We restrict ourselves to the de- 
terministic part and drop the noise part. As we shall 
see, the NMD equation has a structure which is similar 
to the Landau-Lifshitz-Gilbert (LLG) equation, which is 
the fundamental microscopic equation that has been used 
widely to describe the magnetization dynamics of system 
smaller then the exchange length (i.e all excitations are 
uniform in space, hence the system that can be assim- 
ilated to a macrospin). Therefore a joint study gives 
strong insights on the role of transport in the magneti- 
zation dynamics and we start this section with the LLG 
equation. It was understood in the last few years that 
the non-conservation of the spin current could lead to a 
current induced torquoii, magnetic relaxation^ and in- 
terlayer exchange interaction^. In the framework of the 
NMD equation, we will find that the first two physical 
effects can be viewed as different manifestations of one 
single term in the equation. The latter effect however is 
destroyed by the absence of coherence (T > T) in our 
system. 



1. The Landau-Lifshitz-Gilbert equation with spin torque 

The LLG equation consists of a conservative term (the 
magnetization precesses around trajectories of constant 
energy) plus a phenomenological damping term that al- 
lows the system to relax to its equilibrium. Recently, the 
LLG equation has also been used to describe the dynam- 
ics of a thin magnetic layer in presence of spin polarized 
current—. In those systems, the polarized current exerts 
a torque that drives the system out of equilibrium. The 
spin torque can be incorporated into the LLG equation 
through an additional term. 

We start by writing the LLG equation for the mag- 
netization of the nanomagnet in presence of the current 
induced spin torque, as introduced by SlonczewskiS. Here 
we suppose the magnetic grain is connected to two leads, 
one of which being magnetic. See Fig.^for a cartoon of 
the system. We denote by m a unit vector pointing in 
the direction of the nanomagnets magnetization. 



dm 
~dt 



m x 



1 dE 

HSq dm 



dm I 



g(m 



z) z x m 



(40) 

where e the charge of the electron, / the current flowing 
through the nanoparticle and So the total spin of the 
nanoparticle. The magnetic energy E(m) consists of a 
uniaxial anisotropy (along the z-axis) and a Zeeman term 
that couples to an external magnetic field H (which also 
lies along the z direction). 



E{m) = —nSaim • z) 2 — giisSoH • 



(41) 



The terms proportional to the current / corresponds to 
the spin torque as calculated in Ref— . Here again we 



have supposed that the magnetization mi of the polariz- 
ing lead lies along the z-axis. The torque is asymmetric 
with respect to the parallel/ anti-parallel configuration 
(with respect to mi) and with good approximation, the 
torque is modulated by a function g(M ■ z) that takes the 
following form. 



g(M ■ z) = 



1 + qM ■ z 



(42) 



with < q < 1. The system having a cylindrical sym- 
metry around the z axis, Ea. H4Q(l can be rewritten in 
spherical coordinate (#,(/>) (see Fig. and we get, 



9 = —ot(j)sin8 — -^g(cos9) sin9 



> sin 9 = a9 + ?f cos 9 sin 8 + jH sin 9 



(43) 



The angle <fi can be eliminated and we get a closed equa- 
tion for u = cos 9. 



1 



1 



-(1 



2k la 
ajH + a—u + — _ 
h e£>o 1 + qu 



(44) 



where 7 = gpLs/h is the gyromagnetic ratio. Ea. H44|) is 
the effective LLG equation along the z axis. Introducing 
the dimensionless quantities 



H = 



I 



9^bH 

2k 
Hal 

2aneSo 



(45) 
(46) 



we plot its phase diagram in Fig.[S] This diagram was al- 
ready obtained in ReS^ with a different approach. Here, 
we sketch its construction as it will be useful in the study 
of the NMD equation. 

The phase diagram is constructed from the analysis of 
the stability of the various fixed points of the equation 
(for u — R(u), stable fixed points are given by R(u) = 
and d u R(u) < see FigQJ: the right hand side of Ea. l|4^|) 
vanishes for u = ±1 and possibly for up to two other 
fixed points, depending on / and H . The phase diagram 
is then determined by 3 curves, which defines 5 regions 
(See Fig.0 : 



• Stability of u = +1 : 

H+1 + — — >0 
1 + 9 

• Stability of u = — 1: 

H-1+ -!— < 



(47) 



1-9 



• Existence of two other fixed points u* , u* 
(qH - l) 2 - Aql > 



(48) 



(49) 
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FIG. 4: Sketch of the different possibilities for the u curves 
f Eg. 1441 1 leading to the different parts of the phase diagram. 
The full circles indicate the stable fixed points while the empty 
circles the unstable ones. In diagram (SP) and (SP, +1), an 
intermediate point — 1 < u* < 1 is stable corresponding to a 
spin precession state. 




H 



In the following, we will label a region of the parameter 
space by the list of its stable fixed points : for example 
the (—1) phase corresponds to one stable fixed point at 
u = —1, (SP, +1) corresponds to two stable fixed points 
at u — — 1 and an intermediate value —1 < u* < +1. 
We observe that : i) stable and unstable fixed points al- 
ternate between -1 and 1, ii) if u* < u** exists and are 
located between —1 and 1, u* is stable and u** is unstable 
(because of the sign of the right hand side of 14411 1. From 
this remarks and the stability of u = ±1, we find the re- 
gions (+1),(— 1),(— 1, +1) and (SP) of the phase diagram. 
The region delimited by the 2 lines and the parabola 
(shaded on Fig. [5J requires a more detailed analysis : on 
the parabola, we have u* — u** — —(H + l/q)/2 (they 
reach ±1 at the points A and B, where the stability lines 
of -1 and +1 are tangent to the parabola). A straight- 
forward calculation shows that the two others frontiers 
of this region correspond to h* = —1 and u** — 1 re- 
spectively. Hence, in this region, —1 < u*,u** < 1, and 
we have an intermediate stable SP fixed point (and an 
intermediate unstable one) : it is the (SP,+1) region. 
The five possible fixed points configurations are shown 
in Fig. 21 where a sketch of ii is presented. The pres- 
ence of the (SP,+1) region is physically important since 
it removes the possibility of switching the grain's magne- 
tization back and forth without hysteresis. We also note 
that when the current induced torque is symmetric with 
respect to the parallel/anti-parallel configuration (q = 0) 
the region (SP) and (SP,+1) are sent to infinity. 



2. The NMD equation 

In this section, we study the phase diagram of the 
NMD equation l|24l) at zero temperature for various sets 
of parameters, both analytically and from the numerical 
solution of the fixed points equations and analysis of their 
stability (See Fig. [BJ. The phase diagram is obtained at 
zero temperature by analogy with the LLG equation l|44|) 
with the additional complexity that the presence of the 



FIG. 5: Phase diagram of the LLG equation 140H . The 
various phases are denoted by the list of their stable fixed 
points (See text). The intersection points have the follow- 
ing coordinates in the (H, I) plane : A = ^—2 — |, J , 

s =( 2 -i^)^=H^)- 



Coulomb blockade phenomena makes the "current" term 
in Ij24(l more complicated than the one of the LLG equa- 
tion. Indeed this "current" term can be switched on and 
off as u goes through the critical values where the addi- 
tion energies AE a (u) = e — nau — h% are equal to the 
chemical potential of the leads: the conducting channels 
for majority (resp. minority) electrons can open or close 
when u varies. 

The simplest case is the symmetric case, where the 
torque asymmetry q vanishes, this removes the possibil- 
ity of exciting spin precession states (there can only be 
one intermediate fixed point u* and according to the sign 
of the coefficient of u it must be repulsive). In the up- 
per left panel of Fig|SJ we have represented such a case, 
where the two leads are magnetic with same tunneling 
rates, and lie in the anti-parallel configurations. This is 
the most favorable case for current induced switching: 
not only up electrons have a higher probability to enter 
the grain than the down electrons, but in addition, once 
in, they have a lower probability to get out so that the 
imbalance between up and down current is enhanced. 

To study the more interesting case q ^ 0, we will now 
restrict ourselves to the case where the chemical poten- 
tial of the second lead is very low. In the Lower Left 
panel of Figure 10, we display such a case : note the 
presence of two additional regions of stability (SP, +1) 
and (— 1,SP, +1), the latter having no equivalent in the 
LLG equation. To gain analytical understanding, it is 
convenient to rewrite the NMD equation as : 

u = (l-u 2 )(h + u+^^-) (50) 
V l+quj 



14 




- (+1) 



(-U1) 

(-1) 



(-1,SP,+1) 
(SP,+1) 
(-1,SP) 
-(SP) 



^bias 



4 4.5 5 5.5 



^bias 



FIG. 6: (Color Online). Phase diagram of the NMD equation 
as a function of the reduced bias voltage v and magnetic field 
h (i.e. measured respectively in unit of fc/e and 2n/y). The 
tunneling rates are in unit of 8SoaBfi/h and Vbias has been 
translated of the charging energy. The color code corresponds 
to the stable fixed points as indicated on the figure. SP,+1 
and -1 indicate the corresponding stable fixed points u = ±1 
or u* 7^ ±1 for the Spin Precession (SP) state. Upper Left 
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where we rescaled the time by t = (2asn)t, restricted 
ourselves to r\ = 1 and defined : 



h 
2k 



(51a) 
(51b) 



i(M ^8^( (1 ~^ (w) - (1 + ^ (w) ) 



(51c) 



where 6 a (u) = 9(ah + au + v) with a = ±1 and 9 the 
Heaviside function. 9 a (u) is 1 when the channel a is open 
and otherwise, for a given value of u. When all channels 
are open, A(u) takes the value : 



A(u) = 9 



h 



(i-?)rl -(i 



9)ri 



(52) 



We start with the stability analysis of the two fixed 
points u = ±1. It is summarized on Fig. 0: we plotted 
the various regions where up and down channels are open 
and closed, together with the line of stability for (+1) and 
(-1) in thick solid line. The phase diagram is then estab- 




FIG. 7: The upper (lower) thick line corresponds to the limit 
of stability of the u = — 1 (u = +1) fixed point in the tension 
bias-magnetic field plane. The various thin lines correspond 
to the opening of the the spin up/down channel depending on 
the magnetic configuration. 



lished in the following way: first we concentrate of the 
large bias region (the right part of the phase diagram), 
where all channels are open. In that case, A(u) = O so 
we keep A fixed and vary q, in order to use our previous 
results on LLG equation; second we reduce the voltage 
to see when various fixed points are destabilized when 
some channels closes. At high voltage, we find 4 possibil- 
ities, depending on how A compares to the /-coordinate 
of points B and C on Fig [3 : 



e < 



a--?) 2 



: this is the small q region. is lower 



than the /-coordinate of the B point and one can 
not excite (SP) fixed points. The phase diagram is 
similar to the Upper Left panel of Fig. [G] 

( - 1 ~ l? - > < O < 1_q : <d is between the /-coordinate 

q - - q 

of the B and the C point on Fig. Therefore, for 
some magnetic fields, there is a stable (SP) point 
and the phase diagram of the NMD equation is sim- 
ilar to the lower left panel of Fig. El If one starts in 
the (SP) region and lower the bias, the (SP) fixed 
point will be destabilized by the closing of the up 
channel at u = u* , not at u = —1. This happens 
at a lower bias, and leads to the existence of the 
(—1, SP, +1) region. Note that the first inequality 
is equivalent to requiring that h + u + 0/(1 + qu) 
has a negative slope at u = — 1. At the bottom of 
the (— 1,£P,+1) region, the frontier merges with 
the stability line of -1, since u* — —1 on [B,C] on 
Fig. [SJ On the contrary, this region is limited at 
higher magnetic field h, corresponding the frontier 
of region (SP,+l) and (+1) on Fig. |S| : at this 
point, (qh - l) 2 - 4 9 = hence h max = 1 ~ 2 ^ g ; 
moreover, on this frontier u* is not ±1, so the up- 
per limit of the (—1, SP, +1) region does not merge 
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with the stability line of -1. 



1-9" 



< e < 



In this case, O is bigger than 



q _— — q 
the /-coordinate of the C point in on Fig. but 
lower than the A coordinate, and there is two addi- 
tional regions (as a function of the magnetic field) : 
(SP) and (—1,SP). This case is presented on the 
Upper Right panel with a zoom in the Lower Right 
panel of Fig. @ The origin of the (-1,SP,+1) 
region is analogous to the previous case. 

• < 0. In this regime, O is bigger than the 

/-coordinate of the A point in Fig. This regime 
corresponds to rather high polarization of the leads 
and is not illustrated in Fig. [6] 

The (—1, SP, +1) region has no equivalent in the LLG 
equation. In this region, the up spin channel is blocked 
at u — — 1 so that the torque term does not destabilize 
the (—1) phase. For u slightly higher, this channels opens 
up, and stabilizes a SP fixed point (provided the torque 
is assymetric enough). 

The behavior in the regions with multiple fixed points 
is clearly highly hysteretic, and the hysteretic loops in 
magnetic field or bias voltage are very different. As an 
illustration, we present in Fig. 03 two cross-sections of the 
Lower Right panel of Fig.HJl In the left panel, we present 
the current / as a function of Vbi as with H chosen such 
that the system goes through the (— 1, +1), (—1, SP, +1) 
and (SP, +1) phases. We find two different unconnected 
loops: if the system is initially in u = +1, it will stay 
there, no matter the value of Vbi as - On the contrary, if 
the system is initially in a SP state or in u = —1, it will 
switch, hysteretically between the SP and the —1 state. 
In the right panel, Vbias = 5.5 and we plot the current / 
as a function of magnetic field H . In the hysteresis loop, 
the value of u* in the SP states evolves continuously as a 
function of H so that the current no longer takes discrete 
values. When the system is in a SP state, in addition 
to the DC current calculated here, small AC corrections 
are expected. This AC current, though small (~ 1/Sq 
smaller than the main contribution) is strongly picked 
on the precession frequency of the SP state (oc ku) and 
can be used for an experimental proof of the presence 
of the spin precessing state. In this article however, we 
did not calculate this AC current, which would require 
to take into account the off-diagonal part of the density 
matrix. 



3. NMD equation without bias voltage 

In addition to the current induced magnetic reversal 
and spin precession state described in the previous sec- 
tion, the NMD equation also describes the extra relax- 
ation induced by the presence of the leads. It is inter- 
esting to note that the current induced torque and the 
extra relaxation are just two manifestations of a single 



5 
4.5 
4 

<3.5 
3 
2.5 



1 


i 1 i 


i 1 i 


- 


< — > 


- 




< — 

* 






(-1,SW,+ 1) 
i , i 


(SW.+ l) 
i , i 


4.6 


4.8 5 
V,. 


5.2 5.4 



+1 



sw 
-1 




FIG. 8: Current as a function of bias voltage Vbias (left panel) 
and magnetic field H (right panel), measured respectively in 
unit of ft/e and 2k/7. The voltage Vbias is shifted by e = 20. 
The tunneling rates are in unit of 8SootBK/h~ and Vbias has 
been translated of the charging energy. Same system as in the 

12, = 4, r r T ight = s 

8. In the left panel, H = —4.1 and the two 
curves correspond to two ways of preparing the system, see 
text. In the right panel, Vbias = 5.5 
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phenomena while in the scattering formalism, they were 
derived as two different contributions 3,33 . When no bias 
voltage is put across the sample, the NMD equation sim- 
plifies into (omitting the noise term) 



(l-« 2 ) 



as{h + 2ku) + 
sinh((ftu + 



r(i- P 2 ) 

8S (l + qr)u) 
h/2)/kT) 



cosh(e/fcT) + cosh((«;u + h/2)/kT) 



(53) 



In the large temperature limit, this equation takes ex- 
actly the form of the LLG equation (|44|l with an effective 
damping constant a c g given by: 



a e ff = ct B 



16fcTS (l + qrju) 



(54) 



When the leads are magnetic and q =/= 0, their contribu- 
tions to the damping becomes anisotropic: the phase di- 
agram is unaffected but the relaxation rates for u = ±1 
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differ. In the low temperature limit, the lead-induced 
relaxation rate saturates to ± gj~h~ijq^^ depending on 
which (up or down) channels arekil opened, and vanishes 
when both channels are opened. 

The zero bias equation for the current simplifies into, 



±y 1 ~ u2 ( riT 



rjqu 



tpl 

6 



n(A_E T ) — n(AE^) 



(55) 

Ea. (|55|l vanishes at equilibrium where u = ±1 as it 
should. However if the system is initially out of equilib- 
rium, a DC current will flow during the relaxation pro- 
cess. If the system is maintained in a spin precession 
state (with radio frequency magnetic field for instance) 
a rectified DC current will follow. 



VI. CONCLUSION 

The discussion of the experimental relevance of the 
NMD equation and of the possibility of observing some 
of the effects described in this paper was already done 
in Ref|18. Beside its application to ultrasmall nano- 
magnets, Eq. I|24fl is also interesting from the point of 
view of spintronic as it can give valuable insights on var- 
ious effects usually discussed in more macroscopic sys- 
tems. Hence, we conclude this article by coming back 
to a few key points raised in this article and contrast- 
ing them with the usual theoretical treatment done on 
macroscopic systems, (i) There is only one kind of elec- 
trons in our system, responsible for both the magneti- 
zation and transport. The possibility of writing a closed 
equation for the magnetization only comes from the sepa- 
ration of time scales between magnetic and charge degree 
of freedom. This timescale separation itself comes from 
the fact that many (oc So) tunneling events are needed 
to change the magnetization significantly (far from equi- 
librium) or from the smallness of the Clebsch-Gordon 
coefficient (close to equilibrium, see^i). (ii) One usually 
distinguish between two kinds of tunneling events, either 
on a majority or a minority state. However, in the case 
studied in this paper, where the system precesses fast 
around its easy axis, the good quantization axis for the 
spin is the easy axis of the grain not the direction of the 
magnetization. Hence we have to distinguish four kind 
of tunneling events depending on the orbital a (major- 
ity/minority) and the spin on the z axis (up/down), (iii) 
Low energy excitations are easily identified in our sys- 
tem. In particular, we can separate the magnetization 
excitations (excitations of type E) from the spin accu- 
mulation (excitations of type C and D, see FigEJ • The 
latter correspond to electrons going to/from a majority 
state from/to a minority one hence changing the total 
spin S while the former correspond to a change of S z 
only. The NMD equation has been derived in a regime 
where only one one-body state carries the current, hence 
there is no spin accumulation in our system. In this 



case, the spin torque and spin accumulation phenom- 
ena completely decouple, (iv) In contrast, the tunnel- 
ing magneto-resistance is very sensitive to the presence 
of spin-accumulation which gives rise to non-equilibrium 
peaks in spectroscopy experiments, (v) The three (up to 
now) phenomena associated with the (non conservation 
of) the spin current in magnetic systems are the current 
induced spin torque, spin pumping induced relaxation 
and magnetic exchange interaction. It is remarkable that 
the first two are contained by a single term in the NMD 
equation. The exchange interaction however, which relies 
on equilibrium spin current and on spin coherence is not 
captured by the master equation, (vi) The possibility of 
stabilizing spin precession states have a similar origin as 
in nanopillars^: the amplitude of the torque decreases 
when the system becomes out of equilibrium. The phase 
diagram is however more complex as the different chan- 
nels for spin up and down can be opened or closed, (vii) 
The torque term in nanopillars vanishes as the sinus of 
the angle between the polarizing layer and the thin free 
layer magnetizations. In contrast, the torque term in our 
nanomagnet can exist even in the collinear regime and 
the torque vanishes as the sinus of the angle between the 
easy axis and the grain magnetization. Hence, we expect 
drastic changes in our system behaviour when the tem- 
perature becomes lower than the tunneling rates and we 
enter the coherent tunneling regime. 
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APPENDIX A: DIAGONAL APPROXIMATION 

In this appendix, we examine under which conditions 
the solution of an equation of the form 



iY 



X = QX + M n e iu ^X 



(Al) 



n=0 



converges to the solution of the equation 
X = QX 

in the limit u>i — > oo, where X is a complex vector and 
M n , Q are nxn matrices. Eq. l|ll|l for the density matrix 
is a special case of equation <|A1[) . In analogy with the 
scalar case, we define 



exp 



a iu>it 



-Mi 



N 



Y = l[A- 1 X 



i=0 



and arrive at 



Y {Ui} (t) = M {ui} (t)Y {a)i} (t) 
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where Mi Ui \{£) —* Q uniformly in t for {u>i} — > oo. De- 
noting by Yrx, the solution of Y = QY, we then have 
Y{LUi}{t) Yoo{t) for all t. However, it does not im- 
ply that lim^oo Y^t) = lim { limt^oo Ys u . } (t). A 
trivial example where these two limits cannot be ex- 
changed is the (scalar) equation y — —ay. There, 
lim a ^o limt^oo y a (t) = while rimt->oo hm a ^ Va{t) = 1. 
Coming back to the Q matrix, we cannot interchange 
the two limits whenever Q has vanishing eigenvalues. In 



the case we are interested in, the matrix Q describes a 
Markov process and hence has all its eigenvalues strictly 
negative^ but one corresponding to the steady state. 
As the Markov process conserves the trace of the den- 
sity matrix (for all {wj}), we can project the dynamics 
into the corresponding hyperplane where all the eigenval- 
ues of Q are strictly negative and therefore the two lim- 
its can safely be interchanged. We conclude that when 
{u>i} — ► oo, one can neglect the oscillating terms. 
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